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ABSTRACT 

3 Sr 

This Report describes the current best estimate of the Ranger 4 
spacecraft flight path and the way in which it was determined. A com- 
parison with independent information sources confirms the accuracy 
of the orbit based on the Deep Space Instrumentation Facility ( DSIF ) 
tracking of the spacecraft transponder for 10*4 hr. The miss parameter, 
as determined by the transponder tracking, is believed to be within 
30 km of the correct value. This error is well within the bounds expected 
and testifies to the accuracy potential of Earth-based tracking. 


I. INTRODUCTION 


This Report describes the current best estimate of the 
Ranger 4 spacecraft flight path and the way in which it 
was determined. A comparison with independent infor- 
mation sources confirms the accuracy of the orbit based 
on the Deep Space Instrumentation Facility (DSIF) 
tracking of the spacecraft transponder for 10M> hr. The 
miss parameter, as determined by the transponder track- 
ing, is believed to be within 30 km of the correct 
value. This error is well within the bounds expected 
and testifies to the accuracy potential of Earth-based 
tracking. 

Section II describes the DSIF transponder orbit in 
terms of its trajectory parameters near the Earth, in 
trans-lunar flight, and near the Moon. Symbols used and 
definitions of key trajectory quantities are given. 

Section III summarizes the key events in the tracking 
of the Ranger 4 mission and gives a general description 
of the DSIF stations and tracking modes. 


Section IV describes the DSIF transponder orbit 
determination and compares that orbit with information 
obtained by the Atlantic Missile Range (AMR) tracking 
of the C-band transponder in the Agena booster stage. 

While the spacecraft batteries were depleted at IOV 2 
hr after launch, the radio beacon carried within the 
Ranger 4 spacecraft’s payload, the ‘rough landing” cap- 
sule, continued to operate on its own power supply. The 
weak signals emitted from the tumbling capsule were 
tracked at the DSIF stations throughout the mission. 
Valuable data were taken at both Goldstone stations in 
the several hours prior to lunar impact. Both the doppler 
shift records and time of signal loss at the Goldstone sta- 
tions confirm the accuracy of the previously determined 
orbit. The results are presented in this Report, Section V. 

Section VI gives a functional description of the in-flight 
determination of the flight path together with the tech- 
niques used in editing and weighting the tracking data. 
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II. TRAJECTORY DESCRIPTION 


The Ranger 4 trajectory was made up of a pre-injection 
and a post-injection phase. The pre-injection phase con- 
sisted of all powered flight and coast periods from launch 
to injection (burnout of the last booster stage). The post- 
injection phase consisted of the coast period from injec- 
tion to lunar impact. 

The trajectory characteristics of the pre-injection phase 
were obtained from observed flight data in combination 
with nominal flight conditions (Ref. 1). The trajectory 
characteristics during the post-injection phase corre- 
sponded to the DSIF transponder orbit (Section IV-B). 
The miss parameter B was used to measure the miss 
distance for the lunar trajectoiy. The miss parameter B 
is defined in Appendix A. 

A. Pre-injection Phase 

Using the Atlas D/Agena B boosters, the Ranger 4 
spacecraft was launched from the Atlantic Missile Range 
on April 23, 1962 at 20 hr, 50 min, 15 sec (20:50:15) 
Greenwich Mean Time (GMT). The fact that the Ranger 
4 spacecraft impacted the Moon without the aid of a 
midcourse maneuver demonstrated the adequacy of the 
performance obtained from the Atlas and Agena boost 
vehicles. The sequence of events from launch through 
injection is shown in Fig. 1. 

After rising vertically for a short period, the Atlas 
booster rolled to a launch azimuth of 100.4 deg (east of 
north), as determined by the launch time, and performed 
a programmed pitch-down maneuver until the booster 
engines were cut off and jettisoned. During the subse- 
quent Atlas sustainer and vernier stages, adjustments 
in vehicle attitude and engine cutoff times were com- 
manded as required by the ground guidance computer 
to adjust the altitude and velocity at Atlas vernier engine 
cutoff. The protective shroud covering the Ranger 4 
spacecraft was ejected during the Atlas vernier stage. 

After the Atlas/ Agena separation, there was a short 
coast period prior to the first Agena ignition. The Agena 
B /Ranger 4 spacecraft was nearly horizontal throughout 
the first Agena burn. The attitude was maintained by 
horizon scanner instrumentation and gyros within the 
Agena booster. At a preset value of sensed velocity 
increase the Agena engine was cut off. 

The Agena B /Ranger 4 spacecraft continued coasting 
in a circular parking orbit for 254 sec at an altitude of 


185 km and a space-fixed velocity of 7.800 km/sec. The 
parking orbit was terminated by a stored command deter- 
mined by the ground guidance computer and transmitted 
to the Agena during the Atlas vernier stage. 

The second Agena ignition, which terminated the park- 
ing orbit, initiated the final increase in velocity prior to 
injection. During the second Agena burn (as was the case 
for the first Agena burn), the vehicle s horizontal attitude 
and engine cutoff were controlled by the horizon scanner 
instrumentation and the preset value of sensed velocity 
increase, respectively. The second Agena cutoff con- 
cluded all powered flight for the Ranger 4 spacecraft and 
represented the injection time. 

B. Post-injection Phase 

Prior to injection, the Agena /Ranger 4 spacecraft 
traveled in a southeasterly direction over the Atlantic 
Ocean. Injection occurred in the mid-Atlantic Ocean. 
Following injection, the Agena and Ranger 4 separated 
with the spacecraft continuing on its course over the 
South Africa continent. The Agena booster then, in turn, 
performed a programmed yaw maneuver and ignited its 
retro-rocket. The retro-rocket impulse was designed to 
eliminate interference with the spacecraft operation and 
reduce the chance of lunar impact by the Agena booster. 

At injection, the spacecraft was traveling 10.958 km/sec 
in geocentric space-fixed coordinates at a geocentric 
radius of 6,567.8 km. The spacecraft geocentric distance 
(Fig. 2) increased while the space-fixed velocity (Fig. 3) 
was decreasing. This, in effect, reduced the geocentric 
angular rate of the spacecraft in inertial coordinates until 
at 1.5 hr after injection the angular rate of the Earth 
exceeded that of the spacecraft. This caused the Earth 
track of the spacecraft to reverse its direction from in- 
creasing to decreasing Earth longitude (positive easterly). 
The subsequent Earth track of the spacecraft was similar 
to that of the Sun except with a greater change in lati- 
tude. These characteristics are illustrated in Fig. 4, which 
shows the Earth track of the spacecraft from launch to 
25 hr past injection. 

For the Ranger 4 trajectory, injection occurred at 2.89 
deg past perigee of the geocentric conic with a flight time 
from injection to lunar impact of 63.76 hr. During the 
first 40 hr past injection, the spacecraft was for the most 
part under the influence of the Earth’s gravitational field 


2 


JPL TECHNICAL REPORT NO. 32-345 


EVENT 

1. LIFTOFF 

2. ATLAS BOOSTER ENGINE CUTOFF 

3. ATLAS SUSTAINER ENGINE CUTOFF 

4. ATLAS VERNIER ENGINE CUTOFF 

5 SPACECRAFT SHROUD EJECTION 

6. ATLAS AGENA- B SEPARATION 

7. AGENA - B FIRST IGNITION 

8. AGENA -B FIRST CUTOFF 

9. AGENA - B SECOND IGNITION 

10. AGENA -B SECOND CUTOFF 

II. SPACECRAFT SEPARATION 

12. INITIATE AGENA YAW MANEUVER 

13. COMPLETE AGENA YAW MANEUVER 

14. IGNITE AGENA RETRO-ROCKET 


•«C=} 
13 8 14 









8-9 8 10 



Fig. 1 . Sequence of events 


and essentially remained in an elliptical geocentric orbit. 
The trajectory during this period can be described by 
an ellipse having a perigee and apogee distance of 6,564 
and 606,407 km, respectively, an eccentricity of 0.978, 
and an inclination of 29.70 deg to the Earth's equator. 



Fig. 2. Geocentric distance to probe vs time 
from injection 


As the spacecraft approached the Moons gravitational 
field, a transition was made from the Earth to the Moon 
as the predominant force affecting the spacecraft’s flight. 
After this transition, the trajectory can be described by a 
Moon-centered hyperbola. The hyperbola is inclined 13.5 
deg to the lunar equator with the spacecraft approaching 
the Moon’s surface in retrograde motion. Lunar impact 
occurred at 57.53 deg before perigee of the selenocentric 



Fig. 3. Geocentric space-fixed inertial velocity of 
probe vs time from injection 
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Fig. 4. Earth-track Ranger 4 trajectory 
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(Moon-centered) hyperbola. The hyperbola perigee dis- 
tance was 1,271 km (467 km below the Moon's surface) 
and the eccentricity was 1.380. 

A general sketch of the Ranger 4 trajectory from injec- 
tion to lunar impact is shown in Fig. 5. The inertial 
Earth-centered coordinates used are referenced to the 
Earths equator and the vernal equinox direction. In 
addition, the position of the Moon during the Ranger 4 
flight and the direction to the Sun are noted. This sketch 
illustrates how the initially elliptical path of the trajectory 
was altered as the spacecraft encountered the influence 
of the Moon’s gravitational field. 

The probe was in direct sunlight except for a brief 
period following injection. At 3 min past injection, the 
spacecraft entered the Earth’s shadow and emerged 38.8 


min later. The relative position along the trajectory at 
which these events occurred is shown in Fig. 4. The 
angular relations between Earth, Sun, and spacecraft 
from injection to lunar impact are graphically illustrated 
in Fig. 6, 7, and 8. 

The portion of the Ranger 4 trajectory when the 
spacecraft encountered the Moon is shown in Fig. 9. As 
the spacecraft approached the Moon’s surface, it was 
occulted by the Moon 70 sec before lunar impact at an 
altitude of 529 km. The actual impact could not be 
observed from Earth. Lunar impact occurred on April 
26, 1962 at 12:50:00 GMT. The spacecraft impacted the 
Moon’s surface at a velocity of 2.669 km/sec. The impact 
location was 121.3 deg from the Moon-Earth line at a 
selenocentric south latitude and east longitude of 12.0 
and 231.4 deg, respectively. The spacecraft arrived at the 
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Fig. 7. Sun-probe-Moon angle vs time from injection 


DSIF transponder orbit extrapolated back a few seconds 
for comparison at the nominal injection time (Ref. 3). 
Table D-l (Appendix D) is a key to the trajectory print- 
out. Table D-2 contains the definitions of the printed 
quantities. Constants and conversion factors used in all 
Ranger 4 trajectory computations are listed in Table D-3. 


Moons surface in the forenoon of the lunar day. The 
relative position of the impact location to the Sun s termi- 
nator, sub-solar point, and sub-terrestrial point is shown 
in Fig. 9. The variation in the spacecraft's altitude and 
velocity relative to the Moons surface during the last two 
hr prior to impact is shown in Fig. 10 and 11, respectively. 


A detailed study of the Ranger 4 trajectory can be 
made by examination of the trajectory printout presented 
in Appendix B. In this printout the trajectory parameters 
are listed at selected times from the epoch of the DSIF 
transponder orbit to lunar impact. The printouts were 
obtained from the initial conditions corresponding to the 
DSIF transponder orbit using the Space Trajectory Pro- 
gram described in Ref. 2. 


Trajectory printouts provided in Appendix C (a) and 
(b) demonstrate the closeness of the actual conditions to 
nominal flight conditions at injection and lunar impact. 
Printout in Appendix C (a) is the nominal flight trajec- 
tory. Trajectory printout in Appendix C (b) is just the 
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Fig. 10. Selenocentric altitude of probe vs GMT during lunar descent (last 2 hr) 
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Fig. 1 1 . Selenocentric space-fixed velocity of probe vs GMT during lunar descent (last 2 hr) 
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III. THE TRACKING SEQUENCE OF EVENTS 


A. Introduction 

This Section summarizes the key events in the tracking 
of the Ranger 4 and the Agena stage. Part B describes 
the DSIF post-injection tracking of the Ranger 4 trans- 
ponder and the payload “rough landing* capsule beacon. 
Part C summarizes the AMR post parking-orbit tracking 
of the Agena C-band transponder by the Twin Falls Vic- 
tory Ship and the Ascension Island FPS-16 Tracking 
Station. 

To help interpret the results of the analysis of the 
tracking data given in Sections IV and V, Table 1 sum- 
marizes the key events of the launch to lunar impact 
sequence. When comparing the Agena orbit with the 
spacecraft orbit, it is important to note that all DSIF 
tracking after J 2 occurs after event 5 (Fig. 1) whereas 
some Agena C-band transponder tracking occurs under 
the following conditions: 

1. Before I 2f when the Agena rocket motor is thrusting. 

2. Between I 2 and event 5, when the spacecraft and 
Agena are still mechanically attached (the path of 
the combination differs from the final spacecraft 
orbit due to the imparting of about 0.3 m/sec rela- 
tive velocity at mechanical separation). 


Table 1. Review of key event times 


Event 

Date 

GMT* 

Remarks 

Atlas liftoff, L 

April 23 

20:50:15 


Agena stage parking 
orbit injection, 1, 


20:58:33 

l + 498* 

Agena stage translunar 
orbit injection, f* 


21:04:15 

l + 840* 

Reference epoch for orbit 
determination, E 


2104:19 

l 2 + 4 s 

Mechanical separation of 
Agena and spacecraft 


2106:53 

t + 998“, 1, + 158* 

Ignite Agena retromotor c 


21:13:23 

l + 1388* 

Burnout of Agena 
retromotor 


21:13:43 

l + 1408” 

Loss of spacecraft 
transponder due to 
battery depletion 

April 24 

07:22 

l + 10" 32“ 

Loss of capsule beacon 
signal due to occupa- 
tion by Moon 

April 26 

I2:47:46 b 

t + 63 h 57“ 31* 

■Universal time at event. 

Corrected by — 1.25 sec to account for signal travel time to Stations 2, 3. 
c The purpose of the retro- maneuver is to bias the Agena stage off of a 
nominal lunar impact trajectory. The resultant probability of the Agena 
stage impacting is thus significantly lowered. 


3. Between event 5 and 6, when the Agena orbit is 
slightly changed by the mechanical separation 
velocity. 

4. Between event 6 and 7, when the Agena orbit is 
being changed by the retro-rocket thrust. 

5. After event 7, when the Agena orbit has undergone 
significant change from its orbit prior to event 6. 

In using the Agena C-band transponder data it is quite 
important to employ only the data corresponding to the 
desired Agena orbit. 

B. DSIF Tracking of Ranger 4 Transponder 
and Payload Beacon 

1. General Information 

The detailed characteristics of the Deep Space Systems 
employed in the Ranger 4 mission are given in Ref. 4. 
The names and locations of the stations used are sum- 
marized in Table 2. Stations 2, 3, 4, 5 use 85-ft diameter 
antennas whereas Station 1 (The Mobile Tracking Sta- 
tion) has a 10-ft diameter Az-EI mounted antenna. 

Table 3 indicates the nominal periods of visibility of 
the spacecraft to the participating DSIF stations during 
the course of the mission. Note that the view periods are 
labeled according to the day of “rise” and that the “set” 
times are often on the next day. Note that the signals 
may be received from the spacecraft somewhat before 
“rise” and somewhat after “set” times. 

The DSIF tracking modes are defined as follows: 1 

GM-1. Ground receiver tracks the transponder signal 
in the 2-way doppler mode. The transmitting 
station (designated by an integer q) receives 
the return signal and compares it with the cur- 
rent transmitter signal to generate 2-way dop- 
pler. At the present time this doppler is much 
more accurate than that taken in any other 
mode. 

GM-2. Ground receiver tracks the transponder signal in 
the 1-way doppler mode. The spacecraft return 
signal is obtained from a crystal reference in 


Reference 6 plus Section IV of this Report. Times measured from 
“rise” refer to rise time at the receiving station listed. 
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Table 2. Deep space station locations 21 


Station 

Location 

Geodetic 

latitude 

Astronomic 

longitude 

2,3 

Goidstone, California, U.S.A. 

35.4°N 

1 16.8°W 

1,5 

Johannesburg, South Africa 

25.9°S 

27.7° E 

4 

Woomera, Australia 

31.4°S 1 

136.9°E 

» Ref. 5. Locations are approximate. 


Table 3. Nominal view periods at DSIF stations' 1 


Date of 
rise 

Station 

Rise GMT 

Set GMT 

View 

period 

April 23 

1,5 

21:13:45 

09:04 :33 b 

11*51" 

4 

22:03:16 

00:53:58 b 

2*51" 

April 24 

2,3 

08:28:45 

16:58:54 

8* 30" 

4 

13:22:51 

02:26:1 9 b 

13* 03" 


1,5 

21:01:54 

09:38:22 b 

12* 37" 

April 25 

2,3 

08:42:25 

17:31:54 

8* 49" 

4 

13:49:29 

02:36:52 b 

1 2* 47" 


1,5 

21:19:05 

09:45:3 l b 

12* 26" 

April 26 

2,3 

08:44:08 

12:47:46 c 

4* 04" 

3 Based on 

5-deg elevation angle and post-flight transponder-determined I 

1 orbit. Universal time at spacecraft. 



1 b Set occurs 

on the next day after rise. 



c Loss of capsule beacon signal due to occupation by Moon. 



the spacecraft (q = 0). The accuracy of the 
doppler data obtained is limited by unknown 
small changes in the spacecraft crystal fre- 
quency. This doppler is termed 1-way because 
the doppler shift occurs only on the spacecraft- 
to-ground transit rather than in both directions 
as in GM-1. 

GM-3. Ground receiver tracks the transponder signal 
in the 3-way doppler mode. One DSIF station 
is in GM-1 and another station is “listening in” 
on the return signal. The accuracy of doppler 
generated in GM-3 is being determined on the 
Ranger series of flights but is limited primarily 
by variations in the reference frequency of the 
transmitting station. Improvements are antici- 
pated in the stability of the transmitter refer- 
ence oscillators which will make 3-way doppler 
a primary data type in the future. 

GM-4. Ground receiver tracks the capsule beacon sig- 
nal in the 1-way doppler mode. The doppler 
limitations of GM-2 are present and the value 
of angle tracking is degraded because of the 
lower, and varying, signal level of the capsule 
beacon. 


The only doppler data used to determine the Ranger 4 
spacecraft orbit based on transponder data was 2-way 
(GM-1), whereas angular data was used when the sta- 
tions were in either GM-1, GM-2, or GM-3. Angular data 
from Station 1 (Table 2) was rejected because carefully 
calibrated, more accurate, data was available from Sta- 
tion 5. 

2. Transponder Tracking 

Table 4 summarizes the transmitter number q versus 
time during the mission, as well as the acquisition times 
on the first pass. The most critical times are initial acqui- 
sition in GM-2 and initial times in GM-1, in reverse 
order. After that, delays of under 10 min in transferring 
transmitting responsibilities from station-to-station have 
minor effect on the accuracy to which the spacecraft 
orbit can be determined. 

The information in Table 4 is somewhat compressed 
in that the time of transition from q = 0 to q ^ 0 is 
chosen to be the time when the first valid 2-way doppler 
was received at the transmitting station; the transition 
from q 0 to q — 0 is chosen to be the time of the last 
valid 2-way doppler point in that interval. Thus, Table 4 
is more aptly a list of time intervals in which 2-way dop- 
pler was taken. 


Table 4. Transmitter number and acquisition time$ a 


Trans- 

mitter, 

9 

Time interval 

Receiving 
station, i 

Acquisi- 
tion time 
(GMT on 
April 
23, 24 

Remarks 

0 

Launch 

tot, = 21:29:31 

1 

21:13 

Rise - 1" 



5 

21:15 

Rise + 1" 

1 

Mot, - 23:05:21 

1 

21:29:31 

Rise + 16" 



4 

22:23 

Rise + 20" 

0 

Mot, = 23:16:51 



t , - t , = 6" 

5 

1 

t, to f, = 23:35:51 

5 

23:16:51 


0 

U to t 5 = 23:40:11 



s 

^r 

II 

1 

1 

tr, taf.i = 23:40:11 

1 

23:40:11 

f« - h ~ 26" 

0 

t« to t 7 = 00:09:51 



t~ trt —— 4" 

5 

f- tots = 07:20:51 

5 

0009:51 

U - U - 7‘ 11” 

0 

h on 




» Reference 6 plus Section IV of this Report. Times measured from "rise" 
refer to rise time at the receiving station listed. 
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The shifting of the transmitting assignment from Sta- 
tion 1 to 5 and back to 1 and then back to 5 represents 
a successful execution of the preflight plan. After 96 min 
of good 2-way doppler from previously flight-tested Sta- 
tion 1, about 30 min were allowed to try to obtain good 
2-way doppler from Station 5 where no 2-way doppler 
had been available previously, and then the transmitting 
job was handed back to Station 1 while Station 5 data 
quality was evaluated. Within 53 min after the first good 
2-way doppler point was received from Station 5, the 
data quality had been determined to be excellent and 
Station 5 was allowed to re-establish 2-way lock. 

As an example of the interpretation of Table 4, con- 
sider the first time q = 1 appears in the left-hand col- 
umn. From t x to t 2 Station 1 was transmitting; Station 1 
achieved 2-way lock (GM-1) at the time indicated to the 
right of “1” in the receiving station column. Station 4 
achieved lock on the signal in GM-3 at the time indicated 
to the right of “4” in the receiving station column. The 
time from f 2 to t z was spent in transferring the trans- 
mitting assignment to Station 5. 

3. Capsule Beacon Tracking 

At 10 hr and 32 min after launch the spacecraft trans- 
ponder signal was lost due to depletion of the spaee- 
craft’s batteries. For the remainder of the mission all 
DSIF stations tracked the capsule beacon, except for 
short periods of time during which unsuccessful searches 
were made for the transponder signal. Table 5 sum- 
marizes the periods of beacon tracking for the DSIFs. 

C. AMR Tracking 

1. Introduction 

After burnout of the final stage, two AMR stations 
tracked the Agena Stage C-band transponder. The first 
data near final stage cutoff came from the Twin Falls 
Victory (TFV) ship. Shortly after the TFV lost the trans- 
ponder the Ascension Island FPS-16 tracker acquired the 
transponder and tracked through the sequence of events 
described in Section III- A. Figure 12 illustrates the ele- 
vation angles at the two stations for the first 10 min after 
the reference epoch E (injection time + 4 sec). 

2. TFV Tracking 

The TFV ship began tracking during the final bum 
and sent data to Jet Propulsion Laboratory (JPL) cover- 
ing the interval from 21:03:19 GMT (E — 60 sec) to 
21:08:16 GMT (E + 237 sec). The ship was reported ‘on 
station” at 326° 45' east longitude, 13° SS' north latitude 



Fig. 12. AMR stations view periods and data spans 


(astronomic) during the tracking interval. At E the probe’s 
elevation was 12 deg at an azimuth of 280 deg east of 
north, at a range of 745 km. At closest approach to the 
ship the vehicle range and elevation were 332 km and 
42 deg respectively. The azimuth angle to the vehicle 
when it was at the ship’s horizon was 125 deg east. A 
total of 78 data sets (time-azimuth-elevation range), after 
£, was received at JPL. The data received were already 
corrected for ship’s pitch and roll by means of an inertial 
reference and an onboard digital computer. The data 
were sampled every 3 sec and transmitted over Milgo 
165 digital-radio teletype equipment at a rate of 1 sample 
per 6 sec. 

3. Ascension Island Station Tracking 

The Ascension Island tracker (7.9° south latitude, 14.4° 
west longitude) sent real-time data to JPL covering from 
21:14:18 GMT (E + 9 m 59*) to 21:29:12 GMT (E + 24* 
53*). Since these data were taken after the Agena retro- 
rocket maneuver had begun (Table 1), it will not be 


Table 5. Summary of capsule beacon tracking 


Date 

Station 

Acquisition 

GMT 

End of track 
GMT 

April 24 

2 

08:32 

17:03 


3 

09:04 

17:40 


4 

13:52 

01:58 


5 

21:40 

09:25 

April 25 

2 

08:47 

17:30 


4 

14:23 

02:13 


5 

21:40 

09:32 

April 26 

2 

8:46 

12:47 


3 

8:33 

12:47 


f 1 
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discussed here. In non-real time, data covering the 
interval 21:10:42 GMT (E + 23*) to 22:21:00 GMT 

(E + 76 m 41*) were received at JPL. About 90 sec of this 
data was obtained prior to Agena retro-ignition. During 
this 90-sec interval the elevation angle varied between 


43 and 46 deg. The azimuth and range at the first and 
last point where 26° E, 1235 km and 71° E, 1648 km, 
respectively. The data sets were sampled every 6 sec and 
contained the same measurement types as the TFV ship. 
Transmission was over the Milgo 165 equipment. 
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IV. FLIGHT PATH DETERMINATION USING TRANSPONDER TRACKING 


A. Introduction 

The real-time determination of the parking orbit is 
the responsibility of the AMR. Their pre-injection track- 
ing of the Agena vehicle C-band transponder is important 
in establishing the parking orbit and detecting non- 
standard flight conditions. The AMR supplies JPL with 
parking orbit elements and initial acquisition informa- 
tion for transmittal to the DSIF stations and for prelimi- 
nary estimation of the spacecraft injection conditions. 
The primary post-injection tracking of the spacecraft is 
done by the DSIF. 

The pitfalls in utilizing AMR post-injection tracking 
of the Agena transponder were described in Section III- 
A. The primary functions of the AMR post-injection 
tracking coverage are: (1) evaluation of the Agena per- 
formance, (2) detection of non-standard flight path, and 
(3) assistance in improving the convergence of the Orbit 
Determination Program (ODP) when very limited 
amounts of DSIF data are available. 

Our long-range objective is to utilize AMR post- 
injection tracking along with the DSIF tracking data to 
determine the spacecraft orbit. We are currently testing 
the consistency of the two data sources and determining 
the best ways to use this information. In Part B of this 
Section we describe the results of the flight path analysis 
of the Ranger 4 spacecraft as derived from the 10.5 hr 
of DSIF tracking of the spacecraft transponder. Part C 
describes the preliminary results of our investigation of 


the compatibility of the AMR Ascension Island and Twin 
Falls Victory Ship tracking data with the DSIF tracking 
results. The results of the comparison are very encourag- 
ing and suggest the lines along which our procedures 
must be modified to utilize the AMR data in the space- 
craft orbit determination. 

B. Flight Path Determination Using DSIF Tracking 
of the Spacecraft Transponder 

1. Summary of Data Taken 

The complete sequence of tracking events and ground 
tracking modes is described in Section III. Section VI-C 
discusses the estimation method used. Table 6 summa- 
rizes the data points used in the orbit determination. 

Angle tracking data was used whenever the ground 
stations were in GM-1, GM-2, or GM-3 and the “data 
condition” code indicated good data. Only 2-way doppler 
data (GM-1) were used; the reasons were discussed in 
Section III-B1. Table 6 provides a gross picture of the 
performance of the data taking and handling system; 
Column 3 gives the total number of data points taken at 
each station during the life of the spacecraft transponder. 
The editing of the data, described in Section VI-B1, 
allowed the number of points (and percentage of total) 
listed in Column 4 to be used in the final orbit determi- 
nation. Of particular interest is the number and percent- 
age of data sets rejected for bad format or as “blunder 


Table 6. Summary of data used in orbit determination 


Station 

Data typos 

Points received 

Points used 

Bad format rejection 

Blunder points 

Bad data 
condition 

Rejection limits 
on blunder points 

% of received 

% of received 

% of received 

% of received 

% of received 

1 

Mobile tracker 

2-way 

doppler 

881 


39 

2 

137 

3 cps 

100 


4.4 

0.2 

15.6 

4 

Woomera 

Hour angle, 
declination 

87 


15 

2 

35 

0.15 deg 

100 


17.2 

2.3 

40.2 

5 

Johannesburg 

2-way 

doppler 

428 


0.14* 

11 

26* 

1.5 cps 

100 

88.0 

3.3 

2.6 

6.1 

Hour angle, 
declination 

960 

719 

29 

53 

159 

0.15 deg 

100 

74.9 

3.0 

5.5 

16.6 


* Doppler and angles are given on each data message and the entire message is rejected for any format errors or bad condition; thus, the 2-way doppler rejects 
of columns 5 and 7 are a sub-set of the angle rejects listed below them. 
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points.” As discussed in Section VI-C, no attempt is made 
to unscramble data messages containing any format 
errors. “Blunder” points can create significant problems 
in converging on an orbit when very little data is avail- 
able and hence are important in influencing the time 
required to establish our first estimate of the orbit. The 
number and percentage of the points omitted because 
of “bad data condition” are listed in Column 7. When 
the tracking station operators or automatic detectors 
recognized that the data being transmitted would not 
be usable, the data condition codeword reflected these 
situations. This situation occurred when re-tuning the 
ground transmitter to maximize the signal received at 
the spacecraft, when commands were being sent, and 
during the acquisition phase. 

2. Weighting of the Data 

The data weights were assigned in accordance with 
the policy described in Section VI-C. The weighting 
assigned to the data depends upon the sampling interval 
and, for doppler, the counting time and the range to the 
spacecraft. During the flight, the effective noise due to 
variation of the transmitter reference frequency was cal- 
culated from regular recordings of the transmitter fre- 
quency. The noise in the doppler due to this variation 
never became a dominant factor because the oscillator 
performance exceeded specifications and because trans- 
ponder tracking ended prematurely. Table 7 summa- 
rizes the sample and counting intervals and weights used. 

3. Discussion of Residuals 

Once the data points and weights are fixed, the set of 
initial conditions which minimizes the weighted sum of 
the residuals squared is found by an iterative method. 
The physical constants described in Ref. 7 were used in 
the trajectory calculation. Subsequently the influence of 


variation of GM-Earth on the resultant estimator was 
examined (Section IV-B4 below). 

The differences between the vector of all observations 
and the calculated values based on the converged solu- 
tion is called the vector of residuals. Figures 13 through 
29 are the residual plots, by station, vs time for the data 
types used in the final orbit. The detailed analysis of 
the residuals will be published in another report. The 
Station 1 doppler residuals have a parabolic form due 
to the rounding of the data. Note that the oscillations 
in the Station 5 doppler data are due to the regular 
tumbling of the spacecraft which caused a variation in 
the equivalent phase center. The tumbling effect can also 
be seen in the angle data since there was a wide varia- 
tion in return signal strength due to relative nulls in the 
antenna pattern. 

4. Statistics of Data and Orbit Estimates 

a. Tracking data statistics . The root-mean-squared 
noise (RMS) and mean of the residuals for each station 
is given in Table 8 for each data type used. Note that 
the RMS noise and weights of Table 7 differ signifi- 
cantly in most cases. The difference in angle weighting 
is due to the presence of low-frequency mechanical 

deflection of the DSIF antennas. 

\ 

fo. Statistics of orbit estimate ; data noise . The accuracy 
of the orbit obtained depends on the statistics of the 
tracking noise and on the statistics of all error sources 
which influence the orbital estimate. The tracking noise 
statistics are represented by the “equivalent or worse” 
white noise method described in Section VI-C. The 
Ranger 4 ODP does not “solve for” nor directly include 
the effects of deviations in physical constants such as 
GM-Earth and station locations. Table 9 gives the co- 
variance matrix describing the uncertainty in the space- 
fixed Cartesian coordinates at the reference epoch E, 


Table 7. Summary of weights, sample, and count times 


Station 

Data type 

B to B + 80 min 

B + 80 min on 

Sample 

spacing, 

sec 

Count 

time, 

sec 

Weight, 
cps a 
or deg 

Sample 

spacing, 

sec 

Count 

time 

Weight, 
cps* 
or deg 

1 

2-way doppler 

10 

T b 

0.7 

10 

1 

0.7 

4 

Hour angle. 








declination 

60 


0.18 

60 


0.18 

5 

2-way doppler 




60 

50 

0.20 


HA, declination 

10 


0.45 

60 


0.18 


*5 is reference epoch used for orbit determination {Table 6). 1 cps = c/2f = 0.156 m/sec where the transmitting frequency f is 96 X 10 7 cps. 
b Station 1 doppler is counted ±V2 sec centered about the message time. AM other stations time tag the data at the end of the counting interval. 
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Fig. 14. Station 1 residuals (from 22:00 GMT April 23, 1962) 
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GREENWICH MEAN TIME (APRIL 23,1962) 
Fig. 15. Station 1 residuals (from 23:00 GMT April 23, 1962) 
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Fig. 16. Station 1 residuals (from 00:00 GMT April 24, 1962) 
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Fig. 19. Sfation 5 residuals (from 21:00 GMT April 23, 1962) 
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GREENWICH ME AN TIM E (APRIL 23, 1962) 
Fig. 21. Station 5 residuals (from 23:00 GMT April 23, 1962) 
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Fig. 22. Station 5 residuals (from 00:00 GMT April 24, 1962) 


24 






WAY DOPPLER, cps DEC, deg HA, deg 


JPL TECHNICAL REPORT NO. 32-345 


0.200 


0.000 


- 0.200 

0.200 


0.000 


- 0. 200 
2.000 


0.000 


i 

CM 


- 2.000 


mm 
hr 02 


GREENWICH MEAN TIME (APRIL 24, 1962) 
Fig. 24. Station 5 residuals (from 02:00 GMT April 24, 1962) 
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GREENWICH MEAN TIME (APRIL 24, 1962) 

Fig. 27. Station 5 residuals (from 05:00 GMT April 24, 1962) 
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Fig. 28. Station 5 residuals (from 06:00 GMT April 24, 1962) 
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Fig. 29. Station 5 residuals (from 07:00 GMT April 24, 1962) 
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considering only the noise on the tracking data. The 
covariance matrix is given in terms of its normalized 
correlation matrix and standard deviations of the coordi- 
nates. The coordinates at E are given in Section II. The 
lower part of Table 9 gives the corresponding quanti- 
ties in Earth-fixed spherical (defined in Section IV-C 1) 
coordinates. 

The covariance of errors in knowledge of the coordi- 
nates at E may be ‘mapped” to the target region using 
the miss parameter B (Appendix A) and T L , the linearized 
time-of-flight, as measures of target error. T h may be 
considered to represent the flight time to a vertical im- 
pact (the influence of B on the parameter T h is thus 
removed). Table 10 represents the standard deviations 
and correlation matrix in the B*R, B*T, T h system. The 


Table 8. Tracking noise statistics 


Station 

Data types 

No. of 
Points 

RMS 

Mean' 1 



1 

2-way doppler, cps 

703 

0.639 

-0.005 

4 

Hour angle, deg 

35 

0.009 

- 0.001 


Declination, deg 

i 

35 

0.007 

-0.002 

5 

2-way doppler, cps 

377 

0.078 

-0.002 


Hour angle, deg 

719 

0.020 

-0.002 


Declination, deg 

719 

0.012 

-0.002 

“Antenna angle biases were calibrated out. 


axes of the 1-sigma dispersion ellipse are found by evalu- 
ating the eigen-values of the 2X2 covariance matrix of 
B*R, B*T uncertainties. The results are: major semi- 
axis = 22.9 km, minor semi-axis = 14.6 km, and orienta- 
tion of major axis = 149.6 deg CCW from the R axis. 
The standard deviation of actual flight time to the esti- 
mated impact point is 18.8 sec. It was determined by 
using Table 10 data plus the relationships 


— — = 0.656 sec/km and ~ 0 

0B • T 0B • R 

c. Statistics of orbit estimate ; physical constants. The 
only error source which could significantly degrade the 
target accuracies indicated in Section IV-B4b, above, 
appears to be GM-Earth. A systematic investigation uti- 
lizing the technique first suggested in Ref. 4 was carried 
out to determine the sensitivity of our target parameters 
to changes in the assumed GM-Earth as well as to form 
an independent estimate of that quantity from our track- 
ing data. 

Figure 30 shows the variation of the weighted sum of 
residuals squared (on second and third iterations) as a 
function of the fractional variation in GM from its nomi- 
nal value. The minimum is at —0.7 X 10 5 and the stand- 
ard deviation of this estimate is 3.2 X 10~ 5 . A wide range 


Table 9. Statistics of knowledge of injection conditions ignoring physical constant errors 


Space-fixed Cartesian Coordinates at Epoch E 



Correlation coefficients 

Standard deviation 


X 

r 

Z 

X 

r 

Z 

X 0.290 km 

X 

1 

0.384 

-0.106 

-0.378 

0.144 

0.057 

Y 0.384 

Y 


1 

-0.941 

0.701 

0.940 

-0.874 

Z 0.676 

Z 



1 

-0.868 

-0.983 

0.973 

X 0.648 m/sec 

X 


Symmetrical 


1 

0.854 

-0.941 

Y 1.242 

Y 





1 

-0.979 

Z 2.225 

* 






1 



Ear 

th-fixed Spherical Coordinates at Epoch E 


Correlation coefficients 

Standard deviation 








i 


r 

0 

X 

V 

7 

<7 

r 0.135 km 

r 

i 

-0.682 

0.031 

-0.974 

0.724 

0.768 

0 0.0063° 

0 


1 

0.614 

0.518 

-0.784 

-0.971 

X 0.0035° 

X 



1 

-0.177 

-0.054 

-0.466 

v 0.0943 m/sec 

V 




1 

-0.585 

-0.606 

y 0.001 5° 

7 


Symmetrical 


1 

0.867 

(7 0.0136° 

<7 






1 
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Table 10. Statistics of knowledge of target error 
ignoring physical constant errors 


Standard 

deviation 

Correlation coefficients 


B-R 

B-T 

Tl 

B • R 21.0 km 

B-R 

1 

-0.375 

0.697 

B • T 17.0 km 

B-T 


1 

0.273 

Tl 1 8.5 sec 

Tl 

Symmetrical 

1 


of opinions as to the accuracy of our current knowledge 
of GM is available. The most pessimistic figures are 
around 1 X 10 -5 . Thus, while our answer is encourag- 
ingly close to the adopted value, it affords no new infor- 
mation. We shall continue to assume the adopted values 
of Ref. 7 with an uncertainty of 0.5 part in 10\ 

The degradation of the orbit estimate due to a 0.5 X 
10' 5 fractional error in GM-Earth is described in Table 
11. The change in the converged target coordinates esti- 
mate per 10 -5 fractional change in GM is listed as obtained 
from the previously described computer runs. 

Previous studies of the effect of station location errors 
indicate that less than a 15-km target error results from 
station location errors of 10 3 deg in latitude and longi- 
tude and 37 m in altitude. 

We conclude that our estimate of the orbit should be 
accurate about a 22-km 1-sigma circle in the B plane and 
about 33 sec in linearized time-of-flight after allowing 
for uncertainties in the physical constants. Due to the 
favorable correlation between T L and B *T errors (Tables 
10 and 11) the 1-sigma impact time uncertainty is only 
26 sec. 



C. Comparison of AMR and D5IF Tracking Resutts 

1. Introduction 

Section III has described the sequence of events nec- 
essary to understanding the use of AMR tracking data 
as well as a summary of the AMR tracking data avail- 
able for comparison with the DSIF tracking results. Part 
IV-C2, below, summarizes our analysis of the TFV ship 
tracking results. First, we computed the orbital elements 
(see glossary of terms that follows) based on TFV data 
alone. Comparison with the DSIF orbital elements sug- 
gested that the ship’s location required adjustment. The 
TFV orbital elements with adjusted ship’s location 
showed good agreement with the DSIF-determined 
elements. 


We plan to re-evaluate the Ranger 4 flight data using 
a more sophisticated orbit determination program which 
has just been completed. Here, the uncertainties in 
physical constants and station locations will be handled 
in a rigorous fashion in order to obtain a better estimate 
of the orbit and its uncertainties. 


In Part IV-C3, below, we have applied the same 
approach to the analysis of the Ascension Island data. 
The initial disagreement of orbital elements derived solely 
from Ascension data led to a comparison of Ascension 
observation with those calculated on the basis of the 
DSIF-only orbit. The assumption of a 6000-yd bias in 


Table 1 1 . Variation in estimate of impact conditions with changes in GM of the Earth a 


Fractional 

A B • T, 

AB«R, 

A Lot, 

A Long, 


A Impact 

WEStm 

change in GMf 

km 

km 

deg 

deg 


time, sec 

sec 

+ 2 X 10" 

-55.7 

-0.4 

0.39 

— 2.25 


-70 

-108 

+ 6X 10~ 3 

-167.4 

-0.8 

1.26 

-7.13 


-206 

-325 

-6 X 10 3 

+ 169.2 

0.0 

-1.04 

6.32 


+ 222 

+ 326 

-Estimated change in impact conditions per 0.5 X 10~ 5 change in GMs: AB 

• T = — 14.0, AB • B 

= — 0.03, A Impact time = 

— 17.9, ATi, = - 

-27.1. 
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the Ascension range data puts all residuals and orbital 
elements within the range of expected variation. 

We consider that the problems encountered are not 
serious and can be eliminated in the future so that our 
goal of using AMR data in assisting the determination of 
the spacecraft orbit in real time is not far off. Modifi- 
cations in operational and computational procedures at 
JPL are indicated in order to make proper utilization of 
the potential of the AMR tracking data. 

GLOSSARY OF TERMS 

The two sets of orbital parameters used are Earth-fixed 
spherical coordinates and a set of Keplerian orbital ele- 
ments. All elements are referred to true (instantaneous) 
equator and equinox of date. 

Earth-fixed spherical coordinates 

r 0 Earth center to probe distance, km 
<f> 0 geocentric latitude, deg 
X 0 longitude, east, deg 

v 0 speed in Earth-rotating framework, km/sec 

y 0 path angle of velocity, above local (geocentric) 
horizon, deg 

a 0 azimuth angle of velocity, east of north, deg 

Keplerian orbital elements 
a semi-major axis, km 
e eccentricity 
i inclination angle, deg 
a right ascension of ascending node, deg 
to argument of perigee, deg 
«> + v sum of a) and true anomaly at epoch E 

2. Twin Falls Victory Ship Tracking Results 

The TFV data available (Fig. 12) brackets the time of 
mechanical spring-separation of the Agena stage from the 
spacecraft. Since the relative separation velocity is only 
about 0.3 m/sec, the Agena orbit and the post-separation 
Ranger 4 spacecraft orbit were treated as one. As dis- 
cussed in the Introduction, we first determined an orbit 


using TFV data only. The weighting standard deviations 
used were 30 m in range and 0.3 deg in azimuth and 
elevation (by comparison with the RMS residuals it can 
be seen that we assumed correlated errors were present 
in both the range and angle data). Table 12 lists the 
Earth-fixed spherical orbital elements at the reference 
epoch E as well as the RMS error of the residuals. 

Table 12 lists the corresponding orbital elements for 
the DSIF orbit found in Section IV-B. In interpreting 
the differences it is important to note that the ships 
position estimate cannot normally be trusted to better 
than ±2 nautical mi (±0.034 deg). To illustrate the 
effect of ship's location on the orbit estimate, the lati- 
tude and longitude of the ship were varied by 0.1 deg 
in turn, with the following results: 


Change in 
orbit 
estimate 

Change in ship’s location 

0.1 deg latitude 

0.1 deg longitude 

8r 0 , km 

-0.0055 

0 

deg 

0.0993 

0 

SA„, deg 

-0.0030 

0.100 

8u 0 , km/sec 

-0.00002 

0 

Sy„, deg 

0.0006 

0 

8<r 0 , deg 

0.0132 

0 


Utilizing the above information, we made an approxi- 
mate adjustment of the ship's location to better match 
the orbital elements; the latitude was changed by about 
0.10 deg and the longitude by about 0.02 deg. Table 13 
summarizes the results analogous to Table 12 with the 
ship's location adjusted. In columns 5 and 6 we have 
listed the expected l-o- errors in each orbital element due 
to the data noise and a ±2 nautical mi error in ship’s 
location. 

The residuals for the adjusted orbit are shown in Fig. 
31. Note that the RMS of the residuals are essentially 
the same for both the adjusted and unadjusted orbits. 
This is because changes in ship’s location can be so well 
compensated for by errors in the orbit parameters </> 0 , A 0 , 
as illustrated previously. 

The only discrepancy between these two orbits which 
appears significant is the 1.7-m/sec difference in speed. 
We believe that this difference is accounted for by the 
ship’s speed of 5 knots (2.6 m/sec) during the tracking 
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04:20 05:00 05:40 06:20 07:00 07:40 0820 

GREENWICH MEAN TIME 

Fig. 31. TFV (adjusted) residuals 


period. In the future, proper arrangements will be made 
for considering the effect of ship's velocity. 

3. Ascension Transponder Tracking Results 

The Ascension data were placed in the Orbit Determi- 
nation Program and an orbit was found for the Agena 
vehicle (pre-retro) using these weights: azimuths and 
elevation angles 0.02 deg, range 30 m. The resulting orbit 
is given in Table 14 where comparison is made with the 
DSIF orbit. Two coordinate systems are shown. Earth- 
fixed sphericals and Keplerian orbital elements, in order 
to emphasize separate features. Note how conspicuous 
is the excessive semi-major axis and r n associated with 
the Ascension orbit. 


In order to detect possible errors, we calculated the 
Ascension Island observations based on the DSIF orbit. 
The residuals are shown in Fig. 33. It appears that the 
Ascension data had range readings which were 5.5 km 
too high. 

A second Ascension orbit was then computed, the same 
as before except 6900 m were subtracted from the Ascen- 
sion range data. This range adjustment was chosen after 
several tries. The results are shown in Table 15. Residuals 
of this orbit are shown in Fig. 34. The fit to the equations 
of motion is better and the discrepancies between DSIF 
and the range-adjusted orbit are smaller. It should be 
noted that if a timing discrepancy existed between Ascen- 
sion and the DSIF, and this were the only error, then 
comparison of orbital elements would show discrepan- 
cies only in two of the elements, Q and (« -f- v). Such does 
not appear to be the case. 


UJ 

o 


< 

(T 


o> 

Q> 

T3 


3 

1 

N 

< 


O' 

a> 

-o 


% 

UJ 

UJ 



21:10:40 21:11:20 21:12:00 

21:11:00 21:11:40 21:12:20 


The residuals of this orbit are shown in Fig. 32 indi- 
cating errors of a systematic nature. 


GREENWICH MEAN TIME 

Fig. 32. Original Ascension Island orbit residuals 
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21: 

:40 li:00 


10:40 ||:00 


GREENWICH MEAN TIME 

Fig. 34. Revised Ascension Island orbit (range 
bias removed) residuals 


GREENWICH MEAN TIME 

33. Original Ascension Island residuals 
based on DSIF orbit 


Table 12. Ship's orbit based on unadjusted location 
Earth-fixed spherical coordinates at epoch E (E = April 23, 1962 21:04:19 GMT) 


Coordinates 


Ship's orbit (unadjusted) 
DSIF 

Difference 

1-ct ship's orbit errors'* 


Data type 
AZ 
EL 

Range 


to, krn 

4>o, deg 

Xo, deg 

v 0t km /sec 

Yo, deg 

Oo, deg 

6568.7855 

14.4806 

320.2411 

10.5416 

1.6847 

117.3531 

6568.8833 

14.5741 

320.2590 

10.5433 

1.6688 

117.2733 

-0.0978 

-0.0935 

-0.0179 

-0.0017 

0.0159 

0.0798 

0.290 

0.0044 

0.0013 

0.0003 

0.013 

0.034 


No. of points 
78 
78 
78 


RMS 

0.0194 deg 
0.0188 deg 
0.0895 km 


Ship's (estimated) position 
(From Section lll-C) 

13.4948°N latitude (geocentric) 
326.7500°E longitude 


*■ Based on theoretical errors due to the data noise implied by the weighting sigmas used. Errors in ship's location and velocity are not considered. 
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Table 13. Ship's orbit based on adjusted iocatton a 

Earth-fixed spherical coordinates at epoch E (E = April 23, 1962 21:04:19 GMT) 


Coordinates 

to, km 

<t>o, deg 

Xo, deg 

v 0 , km/sec 

Y<m deg 

<To, deg 

Ship's orbit (adjusted) 

6568.7814 

14.5735 

320.2563 


10.54155 

1.6852 

1 17.3408 

DSIF 

6568.8833 

14.5741 

320.2590 


10.54329 

1.6688 

117.2773 

Difference 

-0.0109 

-0.0006 

-0.0027 


-0.00174 

0.0164 

0.0635 

1 -a ship orbit errors' 1 

0.290 

0.044 

0.0013 


0.00027 

0.013 

0.034 

1 -o' ship orbit errors, location errors b 

0.002 

0.035 

0.035 


0.00001 

0.0002 

0.005 

Data type No. of points 

RMS 



Ship's (adjusted) position 


AZ 78 

0.0194 deg 


13.5884°N latitude (geocentric) 


EL 78 

0.0196 deg 


326.7680°E longitude 


Range 78 

0.00897 km 






‘Theoretical noise due to data noise implied by the weighting sigmas used. 
b Theoretical noise due to uncorrelated 2 nautical mi latitude and longitude errors. 


Table 14. Comparison of original Ascension orbit with DSIF orbit 

Earth-fixed spherical coordinates at epoch E IE — April 23, 1962 21:04:19 GMT) 


Coordinates 

r Q , km 

<f>o , deg 

X 0 , deg 

Vo, km /sec 

7o, deg 

(To, deg 

Ascension (original) 

6581.141 

14.6300 

320.2565 

10.5579 

1.5080 

1 17.3326 

DSIF 

6568.883 

14.5740 

320.2590 

10.5433 

1.6688 

117.2773 

Difference 

12.258 

0.0560 

-0.0025 

0.0146 

—0.1608 

0.0553 

1 -a Ascension orbit errors 1 * 

1.32 

0.013 

0.004 

0.015 

0.018 

0.017 

Keplerian orbital elements at epoch E 

Coordinates 

a, km 

e 

if deg 

ft, deg 

<o, deg 

Iw + p), deg 

Ascension (original) 

537,275 

0.987759 

29.7699 

334.9101 

146.5029 

149.4230 

DSIF 

306,500 

0.978588 

29.6988 

334.8797 

146.2298 

149.4764 

Difference 

230,775 

0.009171 

0.071 1 

0.0304 

0.2731 

— 0.0534 

Data type No. of points 

RMS 





AZ 

16 

0.031 deg 




EL 

16 

0.0098 deg 




Range 

16 

0.025 km 




‘Theoretical error expected from data errors consistent with the weights assumed. 
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Table 15. Comparison of adjusted Ascension orbit with DSIF orbit 

Earth-fixed spherical coordinates at epoch E (E = April 23, 1962 21:04:19 GMT) 


Coordinates 

r 0 , km 

<p 0r deg 

Ao, deg 

v ol km /sec 

7 or deg 

Oo, deg 

Ascension (adjusted) 

6569.574 

14.5643 

320.2552 

10.5427 

1.6652 

117.3006 

DSIF 

6 568.883 

14.5740 

320.2590 

10.5433 

1.6688 

117.2773 

Difference 

-0.691 

— 0.0097 

— 0.0038 

—0.0006 

-0.0036 

0.0233 

1-a Ascension orbit errors* 

1.32 

0.013 

0.004 

0.015 

0.018 

0.017 

Keplerian orbital elements at epoch E 

Coordinates 

a, km 

e 

i, deg 

deg 

o), deg 

(a> + v) f deg 

Ascension (adjusted) 

306,277 

0.978567 

29.0134 

334.8381 

146.2739 

149.5135 

DSIF 

306,500 

0.978588 

29.6988 

334.8797 

146.2298 

149.4764 

Difference 

223 

-0.000021 

0.0146 

-0.0417 

0.0442 

0.0372 


Data type No. of points RMS 


AZ 16 0.0097 deg 

EL 16 0.0057 deg 

Range 16 0.0061 km 

“Theoretical error expected from data errors consistent with the weights assumed. 
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V. CONFIRMATION OF THE DSIF TRANSPONDER- BASED ORBIT ACCURACY BY 
TRACKING THE CAPSULE BEACON NEAR LUNAR IMPACT 


A . Introduction 

Our analysis of the doppler-shift data received from 
the capsule beacon just prior to impact has confirmed 
that the orbit determined using the DSIF transponder 
data is consistent with these observations. As indicated 
in Section III-B1, the pass of April 26 at the two Gold- 
stone Deep Space Stations began about 4 hr and 4 min 
before the time the spacecraft was occulted by the 
Moon s leading edge. Figures 35 and 36 show the actual 
doppler values recorded at Stations 2 and 3 during the 
last hour before impact with the Moon. In subsequent 
comparisons, the discrete data points are not shown. 

In Part B of this Section we review the data-taking 
system and the formulae relating the observations with 


the spacecraft orbit, capsule transmitter frequency, and 
ground station bias oscillator frequency. Estimates are 
made of the necessary quantities and the actual observa- 
tions are compared to the values derived from the DSIF 
transponder orbit described in Section IV-B. We show 
that the expected variations in the capsule doppler data 
due to errors in our orbit estimate (Section IV-B4) bracket 
the actual observations, indicating consistency of these 
two information sources. 

In Part C of this Section the records at Station 2 and 3 
which define the time-of-signal loss are shown and dis- 
cussed. Again the deviation in the actual loss time from 
that predicted is well within the expected bounds. Plots 



Fig. 35. Actual recorded data from DSIF 2 
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Fig. 36. Actual recorded data from DSIF 3 


are presented which indicate the sensitivity of the time- f tty ~ capsule crystal frequency at a reference time in 

of-signal loss to deviation in coordinates near impact. cps 


8. Data System 

1. Recordings 


D 0 = bias oscillator drift rate in cps/min 
D t = capsule crystal drift rate in cps/min 


The frequency recorded at the tracking station in GM-4 
(Section III-B1) is f c t and is a function of the beacon 
crystal frequency and the bias oscillator frequency at 
the receiving station. Recordings of f ch were made 
throughout the entire mission as described in Section 
III-B3. The formula used in the orbit determination 
program to calculate f ch at the ith receiving station is 
given in Fig. 37. Thus, 

f ch . = 930.15 X 10* 4- (/ 0 + D a A f . - 0.455 X 10*) 

— (ft + D t At,*) ^1 — + higher order terms^ 

where 

f'bi = calculated capsule beacon frequency in cps 


fi — range rate 


CAPSULE BEACON 
« CAPSULE 


CRYSTAL FREQUENCY 



f cbi 


- RECORDED CAPSULE 
BEACON FREQUENCY 


u — bias oscillator frequency in cps 


Fig. 37. Sketch of f cb system 
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c = velocity of light 

A t? = difference in min between reference time and 
recording time f cb . 

Since the capsule crystal-derived frequency is not pre- 
cisely known, and varies significantly with temperature 
and other factors, a value of the crystal frequency as a 
function of time was obtained by studying residuals at 
different stations. By using the relationship f t = f tQ + 
D ( At iy calculated values of f cb t were generated. These 
values showed close agreement to actual data taken at 
Stations 3 and 5 where bias oscillators were steady. The 
value for drift rate ( D t — +1.54 cps/min) in f t corre- 
sponds to a temperature change of —5.3° F/day. Thus, 
having now estimated f t with data prior to the last two 
hours before impact, we must now include the bias 
oscillator drift into the final calculation. 

Periodically, values of the bias oscillators at the various 
DSIF stations were automatically recorded. Other times 


the values were observed, noted as steady or unsteady, 
and manually recorded. Figure 38 shows the oscillator 
recordings at Stations 2 and 3 during the final pass. Note 
that Station 2 recordings were oscillating quite widely 
and were sparse in the last half hour before impact. 
Station 3’s lack of recordings prior to 11*20™ was due to 
visual observation of a steady oscillator. When it did 
start to drift, the operator switched to automatic record- 
ing and then at ll*55 m to manual recordings. The record- 
ings have been represented by the three lines indicated 
in Fig. 38. Since the Orbit Determination Program (ODP) 
can use but one drift constant D 0 , the solid line repre- 
sents the drift D 0 used for the calculated values. Note 
that the solid line passes through all the automatic 
recorded values and deviates from the manual record- 
ings. Therefore, to bring the manual recordings (dash 
lines) into the evaluation, the actual data cards were re- 
constructed to simulate the difference between recorded 
oscillator frequency and the values used in the ODP 
calculations. 


30,656, 580 [ 


STATION 3 



30,656, 500j- 


STATION 2 


30,656,580 p 


OCCULT ATlON 


STRAY- 


520 h 


' A 

X x f \ 
x * % 

V 



30,656,500 L 


GREENWICH MEAN TIME, hr 


Fig. 38. Bias oscillator frequency vs time (Ranger 4 third pass) 
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Residuals were generated [f ch . (observed) — f c&i (calcu- 
lated) on both the original data and the reconstructed 
data using the converged transponder orbit injection 
conditions and the previously estimated values of f t , D u 
and D 0 (Fig. 39). The Station 2 residuals were in agree- 
ment with Station 3 residuals between 11:30 and 12:00 
GMT where bias oscillator recordings were available. 
Since adequate recordings were not available after this 
time, we shall concentrate on the Station 3 doppler infor- 
mation. In most of the discussion that follows we have 
chosen the solid line (Fig. 38) as the representation of 
the bias oscillator frequency. Figure 40 shows the effect 
of the different assumptions on the calculated values 
fcb V The agreement in both cases is discussed in the 
following paragraphs. The manual recordings are con- 
sidered to be questionable since they were derived from 
a counter in a non-standard patch condition. 


2. Discussion of Results 

In order to interpret the accuracy of the transponder 
determined orbit using either of the two curves in Fig. 40, 
it is necessary to describe the expected variation of the 
observable doppler during the last hour of flight. We 
will use the parameters B-T, B*R, and T L described 
in Section IV-B4 to describe the expected variations. 
Table 16 gives the correlation matrix and standard devia- 


Table 16. Statistics of knowledge of target errors 
including physical constant errors 


Standard 

deviation 

Correlation coefficients 


B • R 

B • T 

Tl 

B-R 21.0 km 

B • R 

1 

— 0.290 

0.391 

B • T 22.0 km 

B • T 


1 

0.646 

Tl 32.8 sec 

Tl 

Symmetrical 

1 


tions associated with our estimate of total accuracy 
(Section IV-B4c). 


We have perturbed these parameters in turn by their 
1-sigma uncertainty. Changes of ±20 km in B*R caused 
no significant change in the doppler curves. Variation of 
T l , the linearized flight time, causes a shift of the time 
axis equal to the negative of T L . The resulting doppler 
curves for ±33 sec variation in T L are shown in Fig. 41. 
Figure 42 depicts the effect of ±20 km variation in B • T 
while holding T L , B*R constant. One additional varia- 
tion was made to determine the change in the doppler 
curve caused by a ±0.2% variation in the GM of the 
Moon (Fig. 43). It should be noted that this effect is 
nearly identical with an error in A T L . 

Figures 41-43 are based upon the extrapolation of 
the automatic recordings of the bias oscillator. Also 
plotted is the curve representing the individual data 
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Fig. 39. Residuals on reconstructed data using manual recordings vs residuals 
on extrapolated automatic recordings (Station 3) 
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I! hr 12 

APRIL 26,1962 GREENWICH MEAN TIME, (hr, min) 

Fig. 40. Calculated beacon data— manual recordings vs automatic recordings of bias oscillator at Station 3 


points taken at Station 3. It appears evident that these 
observations are consistent with their expected variation 
based on the uncertainties in the transponder orbit. No 
more careful comparison can be made at this time until 
our computing programs have added flexibility. 

Referring to Fig. 40, the assumption that the manual 
recordings give the correct bias oscillator frequency leads 
to a moderately different doppler curve. Since the offset 
of the beacon transmitter frequency /> 0 is uncertain to 
±20 cycles, the calculated doppler frequency curves 
have a constant offset uncertainty of ±20 cycles. This 
is due to lack of complete calibrations of the bias oscil- 
lator frequency in the region where the slope was evalu- 
ated. Considering this additional factor, it appears that 
either of the two doppler curves in Fig. 40 is reasonably 
consistent with the predicted values. 

Figures 44, 45, and 46 show some of the residual plots 
for Stations 2, 3, and 5. Figure 44 shows the wide 


oscillations due to the (uncompensated) effect of the 
bias oscillator frequency variations evident in Fig. 38. 
Figure 45 shows the residuals at Station 5 (South Africa) 
during an interval where the bias oscillator was reported 
steady. It can be seen that the slope chosen to represent 
the capsule crystal frequency fits the observations well; 
the bias is due to fixed offset in a reference oscillator. 
Figure 46 shows the Station 3 residuals during the final 
pass. Note that they are stable for the first two hours 
and then start drifting off as indicated by the bias 
oscillator log (Fig. 38). 

C. Verification by Time of Signal Loss 

1. Observational Records 

The primary evidence of occupation of the capsule by 
the Moon is the loss of received signal at the ground sta- 
tion. Various functions related to the received signal are 
recorded by the DSIF on magnetic tape and independ- 
ently on direct-write oscillographs. 
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APRIL 26,1962 GREENWICH MEAN TIME, (hr, min) 

Fig. 41 . Actual data vs perturbations in T L 


Figure 47 is a recording of the receiver functions 
recorded on magnetic tape at DSIF-2 for the last few 
seconds prior to occultation. Figure 48 is a reproduction 
of the DSIF-3 direct-write oscillograph record for the 
same time. It is included to illustrate both types of rec- 
ords from which occultation time was determined. 

In Figure 47, the trace labeled signal strength is the 
one of critical interest. At the time noted by the arrow 
124747, the signal started to decay. The rate of decay is 
characteristic of the 10-sec time constant of the receiver. 

The time associated with the event is determined from 
a binary-coded-decimal (BCD) time code which records 
days, hours, and minutes, and from a 1-pps time code. 
Both the BCD code and the 1-pps code are derived from 
the station secondary standard which is synchronized to 
WWV. In Fig. 48, the BCD code may be seen at the 
top of the trace. In the case of Fig. 47, which is a play- 
back of the magnetic tape, the mechanization of the 
playback recorder precludes the recording of the BCD 
code, so that only the 1-pps code appears at the top of 
the trace. 


In Fig. 48, the channel labeled Acquisition Relay is 
an event channel which marks loss of receiver lock, i.e., 
loss of signal. The time shown on the figure is the time 
of change of state of this relay. It is consistent with the 
time of signal tail-off as shown on Fig. 47. The rate of 
tail-off is lower at Station 3 since the AGC time constant 
is longer, i.e., 300 sec. 

The conclusion is that, neglecting signal time-of -flight, 
the capsule signal was occulted by the Moon at 124747. 
The accuracy to which this time can be determined is 
approximately ±0.3 sec/RMS. 

2. Discussion of Results 

The actual and predicted times of signal loss are: 

Actual 12M7 m 46 s GMT April 26, 1962 
Predicted 12 ft 48™10 s GMT April 26, 1962 

when corrected for light-time. The error of —24 sec is 
consistent with the expected 1-sigma variation of 26 sec 
described in Section IV-B4c. The confirmation given by 
the time the capsule beacon signal was lost significantly 
enhances our confidence in the accuracy of the DSIF 
transponder determined orbit. 
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Fig. 44. Station 2 residuals (from 09:00 GMT April 26, 1962) 
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Fig. 46. Station 3 residuals (from 09:00 GMT April 26, 1962) 
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Fig. 48. Oscillograph recording of receiver functions (Ranger 4 echo DSIF 3) 
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VI. FLIGHT PATH ANALYSIS 


4 . Introduction 

The Flight Path Analysis (FPA) group is the part of 
the Spaceflight Operations team which performs the real- 
time radio guidance calculations as well as the post-flight 
determination of the spacecraft orbit. The functions per- 
formed are depicted in Fig. 49. It should be noted that 
the functions are sometimes carried on simultaneously in 
a single digital computer program. 

B. Operational Description 

1. Data Editing, Analysis, and Evaluation 

Editing, analysis, and evaluation of the tracking data is 
accomplished in several ways. 

a. Teletype (TTY) printed display of incoming data 
are visually scanned in real time to detect any sys- 
tematic errors. 

b. Station reports, both printed and verbal, are ana- 
lyzed to detect any abnormalities. In addition, criti- 
cal information on oscillator drift statistics, frequency 
changes, and changes in transmitter assignment is 
evaluated. 

c. Newly received TTY data is periodically entered, in 
batches, into a large digital computer program called 
the Tracking Data Editing Program (TDEP). The 
TDEP checks the format, data condition code, data 
range, station, and time sequence against the input 
master format and control cards. All data are listed 
along with the reason for rejection of any data point. 
The new data which have not been rejected are 
added to the TDEP’s Master Data Tape which con- 
tains all accepted data. 

d. Once the orbit is reasonably well known, the devia- 
tions of the values of new observations from their 
predicted values (the residuals) are tested to deter- 
mine whether they are within selected rejection 
limits. In this way “blunder points ” are easily de- 
tected before they influence the estimate of the orbit. 

e. The residuals and rejected data points are analyzed 
to determine the validity of the noise models and to 
locate any systematic error source. On the basis of 
the information gained from the evaluation of the 
incoming station reports and tracking data, correc- 
tive action is recommended to the Tracking Director. 


OPERATION AND POLICIES 


2. Orbit Determination 

The tracking data placed on the TDEP s master data 
tape is the basis for forming an Orbit Determination Pro- 
gram (ODP) data tape. Control of the information placed 
on the ODP data tape is exercised through input to the 
TDEP. The ODP and TDEP are linked in such a way 
that the ODP can call the TDEP to add new data to the 
ODP tape. The most important ODP inputs are the edited 
tracking data, the data weights, and rejection limits. 

The policies used in editing the data, and in selecting 
weights and rejection limits, are described in Section 
VI-C. During the flight, new data points are continually 
being added to the ODP tape, weights are revised, and 
residuals from selected converged orbits are plotted and 
printed. The converged ODP output provides an estimate 
of the initial conditions and physical constants (param- 
eters, in general) describing the flight path as well as a 
statistical description of the uncertainties in the parameter 
vector. The estimated covariance of the parameter vector 
is then mapped to other regions useful for interpretation 
of results. Typically, the properties of the “error ellipse” 
in the impact parameter plane (B-plane) are computed 
as well as other quantities useful in considering maneu- 
ver alternatives. 

3. Trajectory Information 

At all times of a typical mission trajectory, information 
is essential to the analysis of spacecraft performance and 
scientific data, to assist in tracking station acquisition, for 
antenna pointing, and for general information. As sug- 
gested in Fig. 1, the basis for forming these trajectory 
estimates varies with the amount of information available 
and is continually updated. 

4. Maneuver Alternatives 

Since the examination of the midcourse and terminal 
maneuver alternatives was not necessary on this flight, 
this function will not be discussed at length. As suggested 
in Fig. 49, the trajectory estimate(s), information on ex- 
pected spacecraft performance and current status, statis- 
tics of current and future knowledge of the flight path, 
spacecraft restraints, and mission objectives dependent on 
the flight path are input into a digital computer program 
which is designed to examine the detailed results of fol- 
lowing the available alternative maneuvers (trajectory 
corrections). Commands necessary to implement any of 
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the various alternate maneuvers are also computed and 
checked. 

C. In-flight Policies 

The JPL Ranger Orbit Determination Program (ODP) 
is designed to find the set of initial conditions at injection 
epoch which causes the weighted sum of squares of the 
residuals (observed minus computed) to be minimized. 
We call our method modi/icd-least-squares (MLS) to call 
attention to the method used in obtaining the weights. 
In the usual least squares (LS) method, the individual 
data points are weighted inversely proportional to their 
expected (or measured) variances in forming the weighted 
sum of the squared residuals. In MLS, the independent 
weighting values are determined by the expected (or 
measured) effective variances 2 . In arriving at the effec- 
tive variance for each data type at each station (vs time), 
consideration is given to the effective correlation width 
of all recognized error sources, the sampling rates, range 
to the spacecraft, counting time, and elevation angle. The 
ODP-calculated covariance matrix of injection errors will 
always give a conservative estimate of the accuracy when 


effective variances (“equivalent-or-worse uncorrelated 
noise”) are used. In editing the data, our policy is that 
it is better to reject a data set with questionable format 
than to attempt the real-time correction of the error. An 
analogous policy is used in weighting the data; there is 
a maximum weight which can be assigned to any data 
point independent of whether it appears that the data 
may be dramatically better in a particular time interval. 
By sacrificing our possibility of extracting the maximum 
possible information during the flight we reduce the sen- 
sitivity to “blunder points” or small “hidden” errors whose 
effect may be very significant. Section IV-B summarizes 
our experience on Ranger 4 in terms of the fraction of the 
received data which was rejected for various reasons. 


2 This approach was first used at JPL by A. R. M. Noton in August, 
1959 in a JPL internal Technical Memorandum 312-522, Effect 
of Correlated Data in Orbit Determination from Radio Tracking 
Data. Further discussion was given by A. R. M. Noton, E. Cut- 
ting, F. Barnes (Ref. 8). T. A. Magness and J. B. McGuire have 
developed mathematical expressions to contrast the performance 
of LS, MLS, and minimum covariance estimators (under JPL 
Contract 950045) in terms of the eigenvalues and eigen-vectors of 
the data noise covariance matrix (Ref. 9). 
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APPENDIX A 

Definition of the miss parameter B 


The miss parameter B is used at the Jet Propulsion 
Laboratory to measure miss distances for lunar and inter- 
planetary trajectories and is described by W. Kizner in 
Ref. 10. B has the desirable feature of being very nearly 
a linear function of changes in injection conditions. 

The osculating conic at closest approach to the target 
body is used in defining B. B is the vector from the 
target’s center of mass perpendicular to the incoming 
asymptote. Let S ; be a unit vector in the direction of the 
incoming asymptote. The orientation of B in the plane 
normal to S 7 is described in terms of two unit vectors R 
and T, normal to S 7 . T is taken parallel to a fixed reference 
plane and R completes a right-handed orthogonal system. 
Figure A-l illustrates the situation. 

Our Ranger 4 work has used the orbital plane of the 
Moon as the reference plane. If W is a unit vector normal 
to the orbital plane (W in direction of R M X V„, where 
Rjf is radius vector to Moon from Earth and \ M is the 
space-fixed velocity of the Moon relative to the Earth’s 
center) then T = S r X W defines our coordinate system. 



So 

OUTGOING 

ASYMPTOTE 


TARGET-CENTERED 

HYPERBOLA 


Fig. A-l. Definition of B*T, B*R system 
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Ranger 4 trajectory printout based on DSIF transponder orbit 

SPACE TRAJECTORIES 

RANGER-4 ORBIT BASED ON TRANSPONDER DATA ONLY 
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APPENDIX C (Cont’d) 

b. Ranger 4 trajectory based on the DSIF transponder orbit printout at 
the nominal injection epoch and at lunar impact only 
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APPENDIX D 


Tables related to trajectory printout 

Table D-l. Ranger 4 trajectory key 


COLUMN 1 
ROW 

2 

3 

4 

5 

6 


1 

GME 

J 

H 

D 

RE 

REM 


GROUP A 2 

G 

A 

B 

C 

OME 

AU 


3 

GMM 

GMS 

GMV 

GMA 

GMB 

GMJ 


INJECTION CONDITIONS TARGET 


JULIAN DATE 

MONTH 

DAY, YEAR 

HR. MIN. SEC. 

GROUP B 4 

GEOCENTRIC 

XO YO 


ZO DXO 

DYO 

DZO 


5 

CARTESIAN 



TO GHA 

GHO 




TIME PAST INJECTION 


JULIAN DATE 

MONTH 

DAY, YEAR 

HR. MIN. SEC. 


GEOCENTRIC 




EQUATORIAL COORDINATES 

6 

X 

Y 

Z 

DX 

DY 

DZ 


7 

R 

DEC 

RA 

V 

PTH 

AZ 


GROUP C 8 

R 

LAT 

LON 

VE 

PTE 

AZE 


9 

XS 

YS 

ZS 

DXS 

DYS 

DZS 


10 

XM 

YM 

ZM 

DXM 

DYM 

DZM 


11 

XT 

YT 

ZT 

DXT 

DYT 

DZT 


12 

RS 

VS 

RM 

VM 

RT 

VT 


13 

GED 

ALT 

LOS 

RAS 

RAM 

LOM 


14 

DUT 

DT 

DR 

SHA 

DES 

DEM 



GEOCENTRIC 



CONIC ORBITAL 

B • T AND B-R 

EQUATORIAL 

COORDINATES 


EPOCH OF PERICENTER PASSAGE 


JULIAN DATE 

MONTH 

DAY, YEAR 

HR. MIN. SEC. 

15 

SMA 

ECC 

INC 

LAN 

APF 

RCA 


16 

VH 

C3 

Cl 

SLR 

APO 

TFP 


GROUP D 17 

TA 

EA 

MA 

DAO 

RAO 

MTA 


18 

WX 

WY 

WZ 

PX 

PY 

PZ 


19 

QX 

QY 

QZ 

RX 

RY 

RZ 


20 

SXO 

SYO 

SZO 

TX 

TY 

TZ 


21 

BX 

BY 

BZ 

MX 

MY 

MZ 


22 

B-T 

B-R 

B 

PER 

OMD 

NOD 



HELIOCENTRIC 




EQUATORIAL COORDINATES 

23 

X 

Y 

Z 

DX 

DY 

DZ 


24 

R 

LAT 

LON 

V 

PTH 

AZ 


25 

XE 

YE 

ZE 

DXE 

DYE 

DZE 


GROUP E 26 

XT 

YT 

ZT 

DXT 

DYT 

DZT 


27 

LTE 

LOE 

LTT 

LOT 

RST 

VST 


28 

EPS 

ESP 

SEP 

EPM 

EMP 

MEP 


29 

MPS 

MSP 

SMP 

SEM 

EMS 

ESM 


30 

EPT 

ETP 

TEP 

TPS 

TSP 

STP 


31 

SET 

STE 

EST 

RPM 

RPT 

SPN 



SELENOCENTRIC 




EQUATORIAL COORDINATES 

32 

X 

Y 

Z 

DX 

DY 

DZ 


33 

R 

DEC 

RA 

V 

PTH 

AZ 


GROUP F 34 

R 

LAT 

LON 

VR 

PTR 

AZR 


35 

LTS 

LNS 

LTE 

LNE 




36 

ALT 

SHA 

ALP 

DR 

DP 

ASD 


37 

HGE 

SVL 

HNG 

SIA 





SELENOCENTRIC 



CONIC ORBITAL 

B-T AND B • R 

EQUATORIAL 

COORDINATES 


EPOCH OF PERICENTER PASSAGE 


JULIAN DATE 

MONTH 

DAY, YEAR 

HR. MIN. SEC. 

38 

SMA 

ECC 

INC 

LAN 

APF 

RCA 


39 

VH 

C3 

Cl 

SLR 

APO 

TFP 


GROUP G 40 

TA 

EA 

MA 

DAI 

RAI 

MTA 


41 

WX 

WY 

WZ 

PX 

PY 

PZ 


42 

QX 

QY 

QZ 

RX 

RY 

RZ 


43 

SXO 

SYO 

SZO 

DAO 

RAO 

TF 


44 

SXI 

SYI 

SZI 

TX 

TY 

TZ 


45 

BX 

BY 

BZ 

MX 

MY 

MZ 


46 

B-T 

B-R 

B 

PER 




GROUP H 47 

XOCTAL 

YOCTAl ZOCTAL 

XOCTAL 

YOCTAL ZOCTAL 




48 


YY MM DDD HH 

TT SS SSS 

SOCTAL 
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Table D-2. Ranger 4 trajectory key definitions 


Group 

Trajectory constant 

Group A 



Row 1 

GME 

Universal gravitational constant times the mass of 
Earth, kmVsec 2 


J 

Coefficient of the second harmonic in the Earth's 
potential function 


H 

Coefficient of the third harmonic in the Earth's 
potential function 


D 

Coefficient of the fourth harmonic in the Earth's 
potential function 


RE 

Earth radius, km 


REM 

Conversion factor for converting lunar ephemerides 
into km, 1 e.r. = 6378.150 km 

Row 2 

G 

Universal constant of gravitation, km 3 /kg sec 2 


A 

B 

C 

Moments of inertia about principal axis for the 
Moon, kg km 2 


OME 

Sidereal rotation rate of the Earth, deg/sec 


AU 

Astronomical unit, km 

Row 3 

GMM 

Universal gravitational constant times the mass of 
Moon, krnVsec 2 


GMS 

Universal gravitational constant times the mass of 
Sun, kmVsec 2 


GMV 

Universal gravitational constant times the mass of 
Venus, km 3 /sec 2 


GMA 

Universal gravitational constant times the mass of 
Mars, kmVsec 2 


GMB 

Universal gravitational constant times the mass of 
Earth-Moon, kmVsec 2 


GMJ 

Universal gravitational constant times the mass of 
Jupiter, kmVsec 2 

Group B 

Injection conditions are vernal equinox cartesian 
coordinates in a geocentric equatorial system. 
The principal direction (X) is the vernal equi- 
nox direction of date and the principal plane 
XY is the equatorial plane of date. Z is along 
the direction of the Earth's spin axis of date. 

Row 4 

XO 

YO 

ZO 

Cartesian components of the probe radius 
vector, km 


DXO 

DYO 

DZO 

Cartesian components of the probe space-fixed 
velocity vector, km/sec 

Row 5 

TO 

Time of injection in seconds past midnight of day 
before launch, sec 


GHA 

HA of Greenwich at injection epoch, deg 


GHO 

HA of Greenwich at midnight of day before 
launch, deg 

Group C 

Inertial position and velocity of the probe. Sun, 
Moon and target body in a geocentric equa- 
torial system. The principal direction (X) is the 
vernal equinox direction of date and the prin- 
cipal plane XY is the equatorial plane of date. 
Z is along the direction of the Earth's spin 
axis of date. Miscellaneous parameters are also 
included. 


Group 

Trajectory constant 

Row 6 X 

Y 
Z 

DX 

DY 

DZ 

Cartesian components of the probe radius 
vector, km 

Cartesian components of the probe space-fixed 
velocity vector, km /sec 

Row 7 R 

DEC 

RA 

V 

PTH 

AZ 

Probe radius distance, km 
Probe declination angle, deg 
Probe right Ascension angle, deg 
Probe space-fixed velocity, km/sec 
Pitch angle of the probe space-fixed velocity 
vector with respect to the local horizontal, deg 
Azimuth angle of the probe space-fixed velocity 
vector measured East of true North, deg 

Row 8“ R 

LAT 

LON 

VE 

PTE 

AZE 

Probe radius distance, km 
Probe geocentric latitude, deg 
Probe East longitude, deg 
Probe Earth-fixed velocity, km/sec 
Pitch angle of the probe Earth-fixed velocity 
vector with respect to the local horizontal, deg 
Azimuth angle of the probe Earth-fixed velocity 
vector measured East of true North, deg 

Row 9 XS 
YS 
ZS 
DXS 
DYS 
DZS 

Cartesian components of the Sun radius vector, km 

Cartesian components of the Sun space-fixed 
velocity vector, km/sec 

Row 10 XM 
YM 
ZM 
DXM 
DYM 
DZM 

Cartesian components of the Moon radius vector, 
km 

Cartesian components of the Moon space-fixed 
velocity vector, km/sec 

Row 1 1 XT 
YT 
ZT 
DXT 
DYT 
DZT 

Cartesian components of the target radius vector, 
km 

Cartesian components of the target space-fixed 
velocity vector, km/sec 

Row 1 2 RS 
VS 
RM 
VM 
RT 
VT 

Sun radius distance, km 

Sun space-fixed velocity, km/sec 

Moon radius distance, km 

Moon space-fixed velocity, km /sec 

Target radius distance, km 

Target space-fixed velocity, km/sec 

Row 13 GED 
ALT 
LOS 

RAS 

RAM 

LOM 

Geodetic latitude of the probe, deg 
Altitude of the probe above the Earth's surface, km 
East longitude of the Sun in coordinate system 
defined in Row 8, deg 
Right ascension of the Sun, deg 
Right ascension of the Moon, deg 
East longitude of the Moon in coordinate system 
defined in Row 8, deg 

“These are Earth-fixed spherical coordinates in a geocentric equatorial sys- 
tem. The principal direction x is directed towards Greenwich and is the 
intersection of the meridian plane of Greenwich with the equatorial plane. 
The principal plane is the Earth's geometrical equatorial plane x y. z is 
along the direction of the Earth's geometrical north direction. 
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Table D-2 (Cont'd) 


Row 14 DUT 


Group D 


Row 15 SMA 


Trajectory constant 

Ephemeris time minus Universal Time, sec 
Adams-Moulton step size, sec 
Radial velocity of probe, km/sec 
Sun shadow parameter, km 
Declination of the Sun, deg 
Declination of the Moon, deg 

Characteristics of the Earth conic in the geocentric 
equatorial system described under Group B 

Semi-major axis, km 
Eccentricity 

Inclination of the orbit plane to the equatorial 
plane, deg 

Longitude of the ascending node, deg 

Argument of pericenter, deg 

Magnitude of the closest approach vector, km 


Row 22 B-T 


Group E 


Trajectory constant 

Projection of the impact parameter B b upon the 
vector T, km 

Projection of the impact parameter B b upon the 
vector R, km 

The magnitude of the impact parameter, 11 km 

Period, min 

Rate of change of argument of perigee, deg/day 

Rate of change of RA of the ascending node, 
deg/day 

Inertial position and velocity of the probe. Sun, 
Moon, and target body in a heliocentric equa- 
torial system. The principal direction X is the 
vernal equinox direction of date and the prin- 
cipal plane XY is the equatorial plane of date. 
Z is along the direction of the Earth's spin 
axis of date. Miscellaneous parameters are also 
included. 


Row 16 VH 


Row 17 TA 


Row 18 WX 
WY 
WZ 


Row 19 QX 
QY 
QZ 


Hyperbolic excess speed, km/sec 
Twice the energy (v/s viva energy integral, 
km 2 /sec 2 ) 

Angular momentum, km J /sec 

Semi-latus rectum, km 

Apogee distance, km 

Time from pericenter passage, sec 

True anomaly, deg 
Eccentric anomaly, deg 
Mean anomaly, deg 

Declination of the outgoing asymptote, 11 deg 
Right ascension of the outgoing asymptote, 1 * deg 
Maximum true anomaly, deg 

Components of a unit vector normal to the conic 

R X V 

W = 

| R X V| 

Components of a unit vector in the direction of 
perigee 

Components of a unit vector perpendicular to the 
perigee direction, vector P, and being in the 
orbit plane Q = WXP 

Components of the unit vector R b 


Row 25 XE 
YE 
ZE 


Row 26 XT 
YT 
ZT 


Cartesian components of the probe radius vector, 
km 

Cartesian components of the probe space-fixed 
velocity vector, km/sec 


Sun probe radius distance, km 
Probe celestial declination, deg 
Probe celestial right ascension, deg 
Probe space-fixed velocity, km/sec 
Pitch angle of the probe space-fixed velocity 
vector with respect to the local horizontal, deg 
Azimuth angle of the probe space-fixed velocity 
vector measured East of true North, deg 


Cartesian components of the Earth radius vector, 
km 

Cartesian components of the Earth-space-fixed 
velocity vector, km/sec 


Cartesian components of the target radius vector, 
km 

Cartesian components of the target space-fixed 
velocity vector, km /sec 


Row 20 SXO 
SYO 
SZO 


Row 21 BX 
BY 


Components of the unit vector S 0 b along the direc- 
tion of the outgoing asymptote 

Components of the unit vector T b 


Components of the impact parameter B, b km 

Components of a unit vector which lies in the 

orbit plane and is normal to the radius vector R. 

m = w x A 


Row 27 LTE Celestial latitude of the Earth, deg 

LOE Celestial longitude of the Earth, deg 
LTT Celestial latitude of the target, deg 

LOT Celestial longitude of the target, deg 

* 

RST Sun-target range, km 

VST Sun-target velocity, km /sec 

Row 28 EPS Earth-probe-Sun angle, deg 

ESP Earth-Sun-probe angle, deg 

SEP Sun-Earth-probe angle, deg 

EPM Earth-probe-Moon angle, deg 

EMP Earth-Moon-probe angle, deg 

MEP Moon-Earth-probe angle, deg 


b See appendix A. 


b See appendix A. 
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Table D-2 (Cont'd) 


Group 

Trajectory constant 

Row 29 MPS 
MSP 
SMP 
SEM 
EMS 
ESM 

Moon-probe-Sun angle, deg 
Moon-Sun-probe angle, deg 
Sun-Moon— probe angle, deg 
Sun-Earth— Moon angle, deg 
Earth— Moon— Sun angle, deg 
Earth— Sun— Moon angle, deg 

Row 30 EPT 
ETP 
TEP 
TPS 
TSP 
STP 

Earth-probe-target angle, deg 
Earth-target-probe angle, deg 
Target-Earfh-probe angle, deg 
Target-probe-Sun angle, deg 
Target-Sun-probe angle, deg 
Sun-target-probe angle, deg 

Row 31 SET 
STE 
EST 
RPM 
RPT 
SPN 

Sun-Earth-target angle, deg 
Sun-farget-Earth angle, deg 
Earth-Sun-target angle, deg 
Moon probe radius distance, km 
Target probe radius distance, km 
Sun-probe-near limb of Earth angle, deg 

Group F 
Row 32, 33 

Inertial position of probe in a selenocentric equa- 
torial system. The principal direction X is the 
vernal equinox direction of date and the prin- 
cipal plane XY is the geocentric equatorial 
plane of date. Z is along the direction of the 
Earth's spin axis of date. 

Row 34, 35, 36 

Selenocentric-fixed spherical coordinates of the 
probe. Sun and Earth in a selenocentric equa- 
torial system. The principal direction X is in the 
direction of the mean Moon-Earth line. The 
principal plane XY is the mean selenocentric 
equatorial plane. Z is along the direction of 
the Moon's mean spin axis. Miscellaneous pa- 
rameters are also included. 

Row 32 X 
Y 
Z 

DX 

DY 

DZ 

Cartesian components of the probe radius vector, 
km 

Cartesian components of the probe velocity vector, 
km/sec 

Row 33 R 

DEC 

RA 

V 

PTH 

AZ 

Probe radius distance, km 
Probe declination angle, deg 
Probe right ascension angle, deg 
Probe space-fixed velocity, km /sec 
Pitch angle of the probe space-fixed velocity vec- 
tor with respect to the local horizontal, deg 
Azimuth angle of the probe space-fixed velocity 
vector measured East of true North, deg 

Row 34 R 

LAT 

LON 

VR 

RTR 

AZR 

Probe radius distance, km 
Probe selenocentric latitude, deg 
Probe selenocentric East longitude, deg 
Probe selenocentric-fixed velocity, km/sec 
Pitch angle of the probe selenocentric-fixed 
velocity vector with respect to the local hori- 
zontal, deg 

Azimuth angle of the probe selenocentric-fixed 
velocity vector measured East of the Moon's 
mean spin axis, deg 


Group 

Trajectory constant 

Row 35 

LTS 

Selenocentric latitude of the Sun, deg 


LNS 

Selenocentric longitude of the Sun, deg 


LTE 

Selenocentric latitude of the Earth, deg 


LNE 

Selenocentric longitude of the Earth, deg 

Row 36 

ALT 

Altitude of the probe above the Moon's surface, km 


SHA 

Sun shadow parameter, km 


ALP 

Illuminated crescent orientation viewing angle, deg 


DR 

First time derivative of the probe radius distance, 
km/sec 


DP 

First time derivative of the probe radius direction, 
deg/sec 


ASD 

Angular semidiameter of Moon as seen from the 
probe, deg 

Row 37 

HGE 

Right ascension of Earth in probe coordinate 
system, 0 deg 


SVL 

Declination of the Moon in probe coordinate 
system, 0 deg 


HNG 

Right ascension of the Moon in probe coordinate 
system, 0 deg 


SIA 

Earth— probe— Moon angle minus ASD, deg 

Group G 


Characteristics of the selenocentric conic in the 
geocentric equatorial system described under 
Group B except centered at the Moon 

Row 38 

SMA 

Semimajor axis, km 


ECC 

Eccentricity 


INC 

Inclination of the orbit plane to the equatorial 
plane, deg 


LAN 

Longitude of the ascending node, deg 


APF 

Argument of percenter, deg 


RCA 

Magnitude of the closest approach vector, km 

Row 39 

VH 

Hyperbolic excess speed, km/sec 


C3 

Twice the energy (Vis viva energy integral, 
kmVsec*) 


Cl 

Angular momentum, kmVsec 


SLR 

Semi-latus rectum, km 


APO 

Apogee distance, km 


TFP 

Time from percenter passage, sec 

Row 40 

TA 

True anomaly, deg 


EA 

Eccentric anomaly, deg 


MA 

Mean anomaly, deg 


DAI 

Declination of the incoming asymptote, 1 * deg 


RAI 

Right ascension of the incoming asymptote, b deg 


MTA 

Maximum true anomaly, deg 

Row 41 

WX 

Components of a unit vector normal to the conic 


WY 

R X V 


WZ 

PX 

PY 

W = 

|R X V| 


Components of a unit vector in the direction of 


PZ 

perigee 


b S«« appendix A. 

f Same coordinate system as defined under B except centered at the probe. 
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Table D-2 (Cont'd) 


Group 

Trajectory constant 

Row 42 QX 
QY 
QZ 

RX 

RY 

RZ 

Components of a unit vector perpendicular to the 
perigee direction, vector P, and being in the 
orbit plane Q = W X P 

Components of the unit vector R h 

Row 43 SXO 
SYO 
SZO 

DAO 

RAO 

TF 

Components of the unit vector So b along the 
direction of the outgoing asymptote 

Declination of the outgoing asymptote, 1 * deg 
Right ascension of the outgoing asymptote, b deg 

Time from injection to epoch of pericenter 
passage, hr 

Row 44 SXi 
SYI 
SZI 

TX 

TY 

TZ 

Components of the unit vector S/ b along the 
direction of the incoming asymptote 

Components of the unit vector T b 

Row 45 BX 
BY 
BZ 

MX 

MY 

MZ 

Components of the impact parameter B, b km 

Components of a unit vector which lies in the 

orbit plane and is normal to the radius vector R 

M = wx ^ 

b See appendix A. 


Group 

Trajectory constant 

Row 46 

B • T 

Projection of the impact parameter B b upon the 
vector T, km 


PQ 

Projection of the impact parameter B b upon the 
vector R, km 


B 

The magnitude of the impact parameter, 1 * km 


PER 

Period, min 

Group H 


Cartesian coordinates and epoch of injection 
conditions in the geocentric equatorial system 
described under Group B. 

Row 47 

XOCTAL 

YOCTAL 

ZOCTAL 

Cartesian components of the probe radius vector 


at injection in octal representation, km 


XOCTAL 

Cartesian components of the probe space-fixed 


YOCTAL 

velocity vector at injection in octal representa- 


ZOCTAL 

tion, km/sec 

Row 48 


Epoch of injection 


YY 

Years past 1900 


MM 

Month 


DDD 

Day of month 


HH 

Hours 


TT 

Min 


SSSSS 

Msec 


SOCTAL 

Sec in octal representation 

The time past midnight Greenwich Meridian Time 
on (DD), month (MM) and year (YY + 1900) at 
which the injection epoch occurs is the time 
determined by the sum of HH, TT, SSSSS, and 
SOCTAL. 


b See appendix A. 
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Table D-3. Ranger 4 trajectory constants and conversion factors 


Constants 


Conversion factors 


Constants 


Conversion factors 


GMsun 
GMvhius 
GMuarth 
GMiSa r th-Moon 
GMm oun 
GMjhrs 
GMj U p j t < r 

Ms un/ Mvt-nus 
Msun / Mtart h 
Mn a r t h / Mm oon 
Msun/MEarth-Moim 
Msun / Mil a r s 
Msun/M J up i t «>r 

Equatorial radius of Earth 
1 AU 

Ellipticity of Earth 
Conversion from feet to meters 
Atmospheric model 
Sidereal rotation rate of Earth 
Universal constant of gravitation 
Speed of light 
Mean Moon radius 


1.32715445 X 10” kmVsec 2 

3.247695 X 10 5 km 3 /sec 2 

3.986032 X 10 5 km 3 /sec 2 

4.03503 X 10 5 km /sec' 

4.900759 X 10 3 km 3 /sec- 

4.297780 X 10' kmVsec 2 

1.267106 X 10 s kmVsec 2 

408645 

332951.3 

81.335 

328908 

3,088,000 

1047.39 

6378.165 km 

1.495990 X 10 s km 

1/298.3 

0.3048 

1959 ARDC 

4.1780742 X 10" 3 deg/sec 
6.671 X 10 20 km 3 /kgsec 2 
2.997 925 X 10 5 km/sec 
1738.09 km 


Moon moments of inertia about 
principal axis 


Lunar and solar ephemerides 


Geometrical Earth model, used in 
locating tracking and launching 
facilities upon the Earth 


A - 0.88746 X 10 29 kg km 2 
B = 0.88764 X 10 29 kg km 2 
C - 0.88801 X 10 29 kg km 2 

The Moon and Sun positions are 
obtained from the joint JPL-STl 
ephemerides. For purposes of 
converting into kilometers, the 
conversion factors are: 

1 AU = 1.495990 X 10 s km 
1 e.r. - 6378.165 km 
Clarke spheroid of 1 866 
a - 6378.2064 km 
b = 6356.5838 km 
e 2 = 0.006768657997291 


Earth potential function: 

GMi; r JRe o H Re 

<P {R, 0) — — |^1 + “jj (1—3 sin" 0 ) + — ——(3 ~ 5 sin* 0 ) (sin 0 ) 

D Ri „ 1 

+ 35 ^ — 30 sin* 0 + 35 sin 4 0 ) 

where 

R = geocentric distance 

0 = geocentric latitude 

1 - 1.62345 X 10 ~ 3 
H = -0.575 X 10 s 
D “ 0.7875 X 10 “ 5 
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